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THE USE OF SOURCE-SINK AND DOUBLET DISTRLBUTIONS 
EXTENDED TO THE SOLUTION OF ARBITRARY BOUNDARY 
VALUE PROBLEMS IN SUPERSONIC FLCW 
By Max. A. Heaslet and Harvard Lomax 

SUMMARY 

A direct analogy is established between the use of soiirce— sink 
and doublet distributions in the solution of arbitrary boundary value 
problems in subsonic wing theory and the corresponding problems in 
supersonic theory. The concept of the "finite part" of an integral 
is introduced and used in the calculation of the improper integrals 
associated with supersonic doublet distributions. The general 
equations developed are shown to Include several previously pub- 
lished results and particular examples are given for the loading on 
rolling and pitching triangular wings with supersonic leading edges. 

INTRODUCTION 

The problem of finding pressure distributions over airfoils of 
arbitrary shape and plan form or of finding airfoils which have 
arbitrary pressure distributions is one of the most fundamental 
problems in aerodynamic theory. At the present time the most 
important and satisfactory approach to problems of this type is 
provided by the methods of so-called thin-alrfoil theory. The 
essential assumptions in this theory are that the perturbation 
velocities induced by the airfoil are small relative to the free— 
stream velocity and that the boundary conditions caji be specified 
in a fixed ref erence -plane . 

Under the assumptions of thin-^irfoil theory the theoretical 
analysis of a problem in wing theory resolves itself into the task 
of determining the solution of a second— order linear partial 
differential equation with prescribed boundary conditions. In the 
case of purely subsonic flow^ Laplace's equation in three dlmeneions 
must be considered, while in purely supersonic flow the differential 
equation which arises is algebraically equivalent to the two- 
dimensional wave equation of mathematical physics. The classical 
solutions of these two equations have been developed along two 
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distinct lines: first ty use of orthogonal functions which can be 
derived In terms of the boundary conditions, and alternatively by 
means of Green’s theorem which In turn utilizes a known particular 
solution of the partial differential equation together with the 
given boundary conditions. 

One partlc^llar solution associated with Laplace's equation and 
subsonic aerodynaailos has been found to be outstanding In Its 
mathematical usefulness and, when Identified with the velocity 
potential, has a physical interpretation which can supply. In direct 
application, added Insist Into the nature of the problem. This 
function Is referred to eis the "fundamental solution" and can be 
developed from the concept of a so— called source. A concomitant 
development to the soTirce potential Is the doublet potential, and 
appropriate distributions of these functions are known to be 
Biifflclent for the solution of all problems in subsonic wing theory. 

The extension of the use of the fundamental solution to problems 
In purely supersonic flow introduces mathematical dlff IciLLtles which 
differ essentially from those encountered at low speeds. Both the 
source and the doublet potentials possess singularities on their 
conical characteristic surfaces or Mach cones and. In the case of 
the doublet, the singularity Is of hl^er order than can be treated 
by elementary mathematical methods. In the historical development 
of the solutions of the wave equation this trouble was circumvented 
by replacing the source potential by other partlciilar solutions of 
the differential equation. As an example, Vol terra (reference l) 
Introduced the integral of the fundamental solution and In that way 
reduced the order of the singularities involved. The ana lytical 
development of Yolterra's theory presents no Inherent difficulties 
(e.g., reference 2) but the phy0lcal slgniflceuioe of the particular 
solution 1s lost, the direct analogy with subsonic theory no longer 
exists, and a certain amount of mathemiatioal inefficiency arises 
since, after using the Integral of the source potential. It Is 
found necessary to resort at the end of the analysis to taking a 
final derivative. 

In this report, following methods Introduced by Hadamard 
(reference 3)* s- general solution to the thin— airfoil problem 
In supersonic theory will be given In tenns of the distribution of 
sources and doublets over the given reference plane. Furtheniiore , 
a discussion of the nature of the boundary values required will be 
given. For properly set problems In wave theory It has been fomid 
necessary to Epee if y, usually, both the required function and Its 
derivative with respect to time along the boundary considered. In 
aerodynamic applications of the wave equation associated with 
lifting surface theory and thickness distributions it will be shown 
that only a knowledge of the imknown function or its normal 
derivative along the boundary Is needed since a relationship between 
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the two functions will he established on the boundarj surface. 

In the theoretical portion of the report a brief presentation 
will be made of the differential eqoiatlons Involved and the two forms 
of the fundamental solution. An outline Is then given of the types 
of boundary value problems encountered and, since the purpose of the 
report Is to extend the concepts of thln-alrfoll theory which are 
\ised In subsonic theory to problems arising In supersonic theory, a 
discussion will be given of the subsonic development as a basis for 
the analogy which exists between the methods of solution corresponding 
to the two regimes of flow. In the dl8C\iBslon of the purely supersonic 
case It will be shown that the Introduction of the concept of 
"finite part" wlU. provide a technique whereby the Improper Integrals 
arising from the use of doublets may be evaluated In a straight- 
forward manner. The applications of the theoretical, developments 
will Include the rederlvatlon of some previously published results 
and will also contain the oal.oulatlon of load distributions for 
rolling and pitching triangular wings with leading edges swept ahead 
of the Mach cone from the vertex of the triangles. 


SIMBOIS 

b span of wing 

c chord of wing 

M free— stream Mach nianber 

n normal to arbitrary surface 

direction cosines of normal, n 
static pressure 
rate of roll about I-axls 
free— stream dynamic pressure 
rate of pitch about T— axis 

fundamental, solution of equation (3) 

[(x-xi)2t<y-yi)2+(z-«x)® 3 
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fundamental solution of equation (4) 

[ (i-xi)®-<y-yi)®-(z-zi)^] 
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artltrary region of Integration 
aurface enclosing region B 

perturbation relooltlea In. direction of I—, T— , and Z-ejces, 
reapeotlrely 

free-atream velocity 

Cartesian coordinates in original space variables 
taroafoirmed system of coordinates 

y 1 m*-iT 

Infinitesimal used In analysis 

aurface along vhlcb stream enters induced field of wing 
oonoimal to arbitrary surface 
direction cosine a of oonormal 

■rorlable representing either acceleration potential, 
velocity potential, or any of the three perturbation 
velocity components 

surface on which boundary condltlona are given 
perturbation velocity potential 

variable representing either acceleration potential, 
velocity potential, or any of the three perturbation 
velocity components 

pressure coefficient 
load coefficient 

/moment about X-aJCla' 
roUlngHnament coefficient ( 


qb X wing area 


hCi/h{Tb/27) 


differential operator ( -—x + -—x + 

\oi oy oz - 
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©rontlal operator 
sign denoting "finite part" 




of integral 


Sutscripts 


subscript denoting value of variable on upper aiarface of 
wing 


subscript denoting value of variable on lower surface of 
wing 


subscript denoting variable of integration 

subscript on r denoting fundamental solution in supersonic 
flow 


Superscript 

superscript denoting value of variable on opposite side 
of T from filed point (x,y,z) 


THEORETICAL DEVELOPMENT 
Linearized Equations and Boundary Conditions 

The linearization of the second— order differential equation for 
compreseible fluid flow is developed under the assumptions of thin^ 
airfoil or small— perturbation .theory. If the velocity vector of the 
free stream is parallel to and in the direction of the positive 
X-axis, the resulting differential equation is expressible in the 
form 


(1-M ) — - + — - + — - = 0 

bZ^ 


( 1 ) 


Tidiere fl represents a velocity potential., acceleration potential, or 
any one of the perturbation velocities viiile M is the constant 
value of the free— stream Mach number. Assuming the plan© of symmetry 
of the airfoil to H© in the XY plane, the boundary conditions 
associated with equation (l) are given for Z=0. Moreover, if 
u, V, and w represent, respectively, the perturbation velocity 
components along the X, Y, and Z-axes, and if the velocity of the 
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free stream is V, the direction cosines of any stream line are 
proportional to the point functions V + u(X, Y, Z), v(X, Y, Z) , 

and w(X, Y, Z) while pressure coefficient — is given hy the 
relation ^ 


T 


2u 


( 2 ) 


Detailed discussions of these resiilts may he found in reference k. 


Introducing the affine transformations 


I = X 

y = -/±(1-M2) Y 

z = *»/±(l-M^) Z 


where the signs \mder the radicals are chosen so that real values 
result, it follows that in the subsonic case (M<l) equation (l) 
reduces to 


dy^ 

while the supersonic case (M >1) yields 

Sx^ Sy^ Sz^ 


(3) 




The fundamental solution associated with equation (3) is 


1 

i = [(x-xi)^ + (y-yi)^ + (z-zi)^] ^ 


or, in terms of the original space variables, 

i = [(x-Xi)^ + 3 ^(y-Yl)^+ p^(z-Zi)^] 


1. 

2 


(5) 


(5a) 


where 


=Cl - M^) 
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When the wave equation is to he considered the fundamental solution 
takes the form 


— = [(i-xi)^ - (7-yi)^ - (z-zi)^] 2 ( 6 ) 

or 

_ 1 

— = [ (X-Xi)2 - p2(Y_Yi)2 - p2(&-Zi)2] 2 

^0 

where 


= (M^ - 1) 

These fundamental solutions represent, respectively, in subsonic 
and supersonic flow the velocity potentials at the point (x,y,z) or 
(X,Y,Z) of unit sources situated at the point (s;i,yi,zx) or 
(Xi,Yi,Zi). The velocity potential of a doublet may be obtained by 
taking a directional derivative of the source potential, the 
direction of the axis of the doublet coinciding with the direction 
along which the derivative is taken. These two functions will be 
seen to be of paramount importance when Green’s theorem is applied 
to the given botindary conditions. 

It remains now to mention the types of boundary conditions which 
appear in problems associated with wing theory. As a convenience to 
the development of the theory the normalized forms (equations (3) 
and (4)) of equation (l) will be used and boundary conditions will be 
assumed known with respect to the x,y,z coordinate system. Retransfoi^ 
matlon totheX,Y,Z system of ajces can be made quite simply wherever 
needed in application. In order to define the boundary conditions, 
two subscripts will be introduced: the first, u, denotes the value 

of the required function on the upper surface, that is the limit of 
the function as z approaches zero from the positlye direction; the 
second, 1, denotes the value on the lower surface, that is, the 
limit of the required function as z approaches zero from the negative 
direction. 

Using these definitions the three boundary value problems of 
principal Interest can be defined as follows: 

1. Symmetrical nonlifting airfoil.— In this case w^ = w^ = 0 

over all of the xy— plane except for the region occupied by the 
airfoil where = -2wj = Aw = f(x,y) the function being determined 

by the geometry of the wing. Over all of the xy— plane, Au = 0. 

2. Lifting plate with specified loading.— It is given that 
Au-= Uu — "uj = 0 over the xy— plane except for the region occupied 
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airfoil where ^u. — f(x,y), ■bh© function heing detenained by 
the specified loading. Moreover, Aw = 0 everywhere, 

3. Lifting plate with specified camber, twist, and angle of 
incidence,— Over the xy— plane Aw = 0 everywhere. And, except for 
the region occupied by the airfoil, Au = 0. Over the region occupied 
by the airfoil w = f(x,y) where f(x,y) is determined by the given 
camber, twist, and angle of incidence. 

It should be pointed out that the problems considered here differ 
from the usual type of boundary value problem encountered. In the 
so-called Dirichlet or Feumann problems, which arise in connection 
with Laplace's equation, the value of the normal derivative of the 
function or of Idle function itself is specified along the boundary 
while the Cauchy problem for second-order partial differential 
equations Involves the knowledge of both the function and a derivative. 
Except for one case in the aerodynamic problems listed above, no 
absolute vaJ-ues are given but rather the Jump in the value of the 
function along the boimdary is prescribed. 


Boundary Talue Problems in Purely Subsonic Flow 

Since the purpose of this report is to extend the concepts of 
thin— airfoil theory which are used in subsonic theory to problems 
arising in supersonic theory, some discussion of the former will be 
given to provide lucidity as well as to furnish a basis for the 
analogy which will be shown to exist between the methods of solution 
arising in the two regimes of flow. 

The method whereby the solutions of the given problems can be 
effected is provided by Green's theorem which relates a volume 
Integral over a region E to a surface integral over the surface 
S enclosing E. If a, fi are any two functions which, together 
with their first and second derivatives, are finite and single— valued 
throu^out E, then for the subsonic case 


( 7 ) 


introduced and the directional derivatives on the left side are 
taken along the normal n, drawn Inward, to the surface S. 
Identifying now the function a with the fundamental solution — 

and specifying that E satisfies Laplace's equation, equation (7) 
simplifies to give 


f. 


8c 

0 G — 

on 8n- 


dS = 


[av^G — G V^®]dE 


^2 n2 

where the Laplaclan operator, 
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( 8 ) 


iirtiere 


p =[ (x-ix)^ + (y-7i)^ + (z-zi)®r 2 


The variables of integration in the equation are ii, jxt ^if while 
I, y, z are the coordinates of a point P either inside or outside of 
the region of integration. 

If the point P is assumed to he Inside the region of Integration, 

it is evident that the function ^ becomea infinite at P, and it is 

r 

necessary to exclude this point from the region if formula (8) is to 
apply. Describing a spherical surface S vlth radius e about the 
point P, and cansiderlng the integral over the two surfaces D and S 
which enclose the region, it can be shown that in the limit as 
€ — > 0 equation (8) becomes 


n (i,y,z) 



( 9 ) 


The physical significance of this last relation follows immediately: 


the tern — 
r 


represents a fluid source and the term 


8(l/r) 

8n 

represents a doublet with its axis lying along the normal to S, both 
source and doublet being situated at the surface point Xi, yi, z^. 

The value of the function fl at the point x, y, z is therefore given 
as an Integral of a source and doublet distribution, the strengths of 
the two being determined directly from the respective boundary values 


of D and 




Equation (9) 


expresses the value of fl in terms of the surface 


values of and 

knowledge of both 
of fl . As can be 
which relates the 
case where P lies 


^ but this relation does not imply that a 
on 

these variables is necessary for the determination 
shown easily, another condition may be established 
two surface vaJ-ues. Applying equation (8) to the 
outside the region of Integration, it follows that 


the integral, is equal to zero and that fl and — on the surface, 
- 8n 

are fimctionally dependent. 
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Sufficient Information is now at hand to provide a solution for 
the thin-alrf oil-houndary value problems. Consider the region R 
bounded by the sy— plane and a hemispherical dome of infinite radius 
lying above this plane. For all problems to which the results will 
be applied, the value of 0 may be assumed equal to zero at infinity. 
The contribution of the surface integral over the hemisphere is thus 
zero and, from equation (9) 




djci dyi 


\rtiere the integration extends over the entire plane. The dlj^ectlonal 
derivatives are necessarily in the direction of the positive z-axls 
and subscripts are again introduced to denote conditions on the upper 
side of the plane. Keeping P fixed and integrating over the lower 
side of the ly-plane, it follows that 


0 




dxi dyi 


where the negative direction of the normal may be ignored since the 
Integral is equal to zero. Subtracting these two equations gives 
the e 35 )re 8 slon 


a (x,y,z) = - 1 


dxi dyi (10) 


the integral extending now only over the area t for which the 
Integrand does not vanish. Equation ( 10 ) is the basic eqmtion from 
which all solutions in subsonic wing theory will be developed. It 
should be pointed out that the derivation proceeded from the 
assun^tlon that the point (x,y,z) lay above the xy— plane. When 
(x,y,z) lies below the xy— plane, however, the derlTatlcai can be 
carried through in exactly the same manner. Such a development 
reveals that equation (10) is general and that no restriction need 
be iiiqposed on the position of (x, y, z) relative to the reference 
plane. 

As a partictilar application of equation ( 10 ) consider a thin 
Bymmetrlcal airfoil at zero angle of attack and set ® where 
is the perturbation velocity potential. Conditions of symmetry 

demand that (1. = il>^ = il >7 = j), while = w„ and 


= WJ 
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Thus, If — vj = iiar 



( 11 ) 


and the Telocity potential Is given by a distribution of source 
potentials. This distribution can be immediately related to the 
slope of the basic section by means of the equation 


!I5ie symmetric airfoil can aJ.so be treated by replacing fl by 
the perturbatlott velocity v and in the case of the thin lifting 
surface with given loading the function ft can be set equal to u, 
Snploying, respectively, condltiona of symmetry and Irrotationality, 


Boundary Yalue Problems in Pxirely Supersonic Flow 

Applications of Green* a Theorem .— The problem to be discussed 
at this point is the extent to which an analogue to equation (10) can 
be developed for supersonic flow fields. The first step in the 
presentation is, once more, the introduction of Green's theorem for 
equation (l) after it has been modified to the form given by equa- 
tion (h) . Employing the operator 



oz oz 

equation (10) becomes 




( 12 ) 




Green's theorem now becomes 
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where is the so-called conormal to the surface S and has 
direction cosines equal to t>. , n_, such that 

'^2 = "2^ ^3 = “3 (14) 

where ni, ris, ns are the direction cosines of the normal to the 
surface S (fig. l) . (The conormal at any point Xi, yi, zx of a 
surface is the mirror image in the plane x = Xi of the normal 
through the same point.) If a and n are perfectly arbitrary 
functions, aside from satisfying the usual conditions of single— 
valuedness, etc., equation (I3) represents an identity and this, fact 
■will be useful at a later time. For immediate purposes, however, cr 
and will be chosen as solutions of the differential equation under 
consideration so that 

= 0 


and, consequently. 



The use of equations (I3) and (I5) depends upon an understanding 
of the physical nature of supersonic flow fields. The essential 
feature of such flow is the presence of Mach cones which correspond 
to the characteristic cones arising in the mathematical study of the 
wave equation. In accordance with these concepts a dlsturhance in 
the flow field can affect the flow only within its aftercone, that 
is, the cone with vertex at the point of disturbance and with axis 
extending in the direction of the undisturbed stream velocity 
vector; conversely, a point in the flow field can be affected only 
by disturbances which emanate from points within its forecone. 

When the disturbances are generated by a wing it is, moreover, 
necessary to speak more specifically about the nature of the leading 
edge of the wing. For all cases considered here the assumption will 
be made that the plan form is a polygon, that is, is composed of 
straight line segnents. If the wing is swept ahead of the foremost 
Mach cone, the cones arising at the leading edge will have as 
envelope a wedge-shaped surface passing through and extending back 
from the leading edge, while if the wing is swept back of the 
foremost Mach cone this cone will be the surface along which the 
air first experiences perturbatiens or disturbances. Thus, a 
point P with coordinates x, y, z is affected by all disturbances 
lying within its forecone F and at the same time behind the forward 
surface \, the nature of the latter surface being dictated by the 
leading edge. In figures 2(a) and 2(b) these surfaces, along with 
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the disturbance plane t, are indicated for two different wing plan 
forms. In the applications of equation (I 3 ) the volume integraJ. is 
limited to the portion of space common to the surfaces F, A, and t 
and the surface integral Involves a discussion of conditions on these 
surfaces. 

Up to this point the analogy between the subsonic and supersonic 
cases^ insofar as the use of Green* s theorem is concerned, is quite 
apparent. The principal difference which occurred was brought about by 
the use of the true normal in the subsonic field together with the fact 
that the xy-plane was covered by a hemispherical dome of infinite radius 
whereas, in the supersonic field, the concept of the conormal was 
Introduced and the volimie to be considered was that enclosed within 
O' Finite region. In continuing the analogy, however, far more 
fomidable obstacles arise. To begin with, the discussion of a and 
n over the surface in the subsonic case was relatively simple. Thus, 
with no limitations of generality n could be assumed zero at 
infinity and An was specified completely in the xy— plane. But in 
the supersonic case, althou^ An can be assumed known in the xy-plane 
and, as will be seen later, n may be evaluated on the forward boundary 
of the region, nothing is known of n on the forecone F. Hence a 
must be chosen properly so that the knowledge of n is unnecessary on 
F. The most obvious choice of o would be a particular solution of 
equation (4) which woiUd make a = 0 on F and this is in fact the 
choice used by Tolterra (reference l) and applied to aerodynamic 
problems in reference 2. However, if the analogy is to be maintained 
the choice of a is not arbitrary but must be the three-dimensional 
supersonic source corresponding to the fundamental solution 

- in subsonic theory. But such a solution, — = 

^ rc 

[ — ( 7 —yi)^ — (z— Zi)^] 2 becomes Infinite along the 

forecone F which has the equation (x— Xi)^ — (y— yi)^ — (z— Zi)^ = 0. 

It is Just this difficulty which apparently invalidates any further- 
ance of the analogy and the prediction in advance of an aerodynamic 
shape from a distribution of sources and doublets in supersonic flow. 
However, it is also precisely this difficulty which is overcome by 
Hadamard* s general methods. 


Extension of analogy by Hadamard*s Method.— The full develop- 
ment of Hadamard* s methods cannot be given here, but a rough sketch 
of his reasoning will perhaps be useful. The basis of his arguments 
stems from equation (13)* First it is admitted that the right-band 
side of equation (13 ) will tend to infinity as the surface S 
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approaches T so that ^ is not a regular solution to I = 0 

on r . However, as has been mentioned, equation (I3) still must hold 
whether or not a or fl satisfy the wave equation and thus it still 
provides an equality. Hence, 'if the surface integral tends to 
infinity so also must the volume integral. Further, equation (I3) 
implies that these infinite portions Just cancel since the difference 
of the two integrations must always give zero. To deal with such a 
problem quantitatively by the usual mathematical techniques would 
require the study of a limiting process for each new boundary value 
problem, Hadamard’s contribution was the introduction and Justifica- 
tion of a concept \riiich removed the necessity for studying the 
infinite portions Involved, This concept is best presented by 

means of a new notation, thus the sign | is used and is 

to be read "the finite part of," 

Using this concept it 1 s possible to show that if a were set 

eqtoal to — , then equation (13) could be written 
^c 



ff 


)-i- 

dS 


So \r(/ 

rc 



(16) 


BO that the "finite parts" of each side of the eq\iation would be 
equal. Such a notation would, of course, in general be meaningless 
since in discarding arbitrarily a part which tended to infinity it 
would be possible, by proper combinations, to obtain as a remainder 
any finite value. The fact is, however, that the integrals Involved 
in equation (I6) tend to infinity only at a limit of the Integration 
and this limit always involves the forecone F. It was consequently 
possible to devise a manipulative technique to handle equation (I6) 
without considering the singularities individually. It might be 
mentioned, without stressing the correspondence, that a treatment of 
improper Integrals is also employed in the use of Cauchy's principal 
value. In subsonic thin-alrfoll theory and lifting-line theory 
integrals of the latter type are well known in the form 
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T r° f(xo)dJCo ^ ^ 

li = / o< z<c 

Jo X-Xq 

Ix certainly tends to Infinity as Xq approaches x but the use 
of Cauchy’s principal value allows the very large values of the 
integrand obtained when Xq is on either side of x to Just cancel 
in such a way that Ix is finite and unique. So again the 
integral 


I 2 



A(x)dx 

(Xq-x) 


is finite and vinique and given by Hadamard in the form 
T- A(x)-A(xo) 2A(xo) 

It is actiially possible to generalize the idea of "finite part’^ to the 
case when the exponent in the denominator is of the form J + 2 where 
J is a positive integer but such a generalization is not needed for 
aerodynamic applications and will therefore be omitted. 

In actual calculation, the evalmtion of the integral Ig can 
be shortened considerably. Thus, if the indefinite Integral of 


P A(x)dx 


is written in the form f(x)+C then 


I 2 = 




A(xx)-rA(xp) 


dxx 


2 A(xq) 

(Xp^)VsJ 



F(x)+C-^’(a)-€- 


2 A(xq) 

(xo-x)^^®- 


It follows that if C is chosen so that 

2 A(xq) 


C = 


11m 




— F(x) 
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then the expression for la may be written I2 = -[F(a)+C] 

When C is chosen in this maimer, the notation for the calculation 
may thus be modified to the form 


A(x)dx p * A(x)dx 

where the asterisk indicates that the upper limit is not substituted 
into the Indefinite integral^ F(x) + C* 

The technique for the calculation of the finite part has 
therefore been reduced to three slaple steps: first, the indefinite 
integral F(i) + C is determined, second, the constant C in the 
Indefinite Integral la evaluated by means of a limiting process, 
third, the lower limit of the Integral is substituted into the 
indefinite integral and a minus sign prefixed. As an example, 
consider the integral 


-'a (xo2-x2)3/2 


'Jl (Xo2-5) 3/2 


In this case F(x) + C 


— — II ■ — -f 0 and 



- { 2 io) 


= 0 


so that, finally. 



xdx 

(Xq^xS) 3/2 


-[F(a) + C] » 

(xo2-a2)i/2 


With the aid of this artifice the analogy between the subsonic 
and supersonic cases can be continued with relative ease. Thus, in 
equation (I6) the ri^t-hand member is zero provided we exclude the 
point P from the volume of integration. This can be done most easily 
by limiting the integration to the xi = constant plane, a distance e 
upstream from P. The portion of this plane intersected by the cone, 
and thus the section over which the integration must be carried, will 
be denoted by z (fig. 3). 
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Aa drawn, figure 3 " shc3ws a cross section in a yi = constant 
plane for the special case when P is located directly behind and 
above the foremost distvirbance . Applying equation (16) to the 
regions above and below the dist\irbance surface x (plane of the 
airfoil) yields the two equations 


If 


A+r+x'- 




1 " 
rc Bu - 


dS 



1 . ^ 


dS 


(17) 


and 


ff 


x+r+x 



_i_ ^ 
rc 


dfl = 0 


(18) 


where the prime indicates the surface value of ft on the opposite 
side of T from P. 

The integration over Z can he computed for e very small. 

For convenience^ consider P to he the origin; then it follows that 
since the cononnal is in the Xi direction and the area element can 

he written y dy d 6 where 6 = arc tan ~ and y = "tii© 

71 

ri^t aide of equation (I7) will give 
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n 



i_ 

rc ^ 


dS 


11m 

e — > 0 




^ -c y dr 
[e2 - y2]3/2 


^fl(x,yyz) 


hxi J 

0 Jo /^2 — r2 


- llm s27rfl(i,y,z)| —5 I_Z_ ^ 27 t 

e _> o'- ^ [e - y ' 


^n(x,y,z) /^e 7 


-/ 

'-'n 


^ /eS _ 7^ 


2TTil(xjy,z) 


Hence the value of n at the point P, 0.{x,j,z) can he determined 
from equation (I7) with the restriction implied hy equation (I8). 
Further, since only the "finite part" is considered, the integration 
over r yields zero and the two equations comhlne to give 


ft — — — 

2rt 



1 '6 f 1 ,, 1 

— Vr r — )~'T~\ — )(ft^') dxi dyi 

r« ' ^z ^z ' ^z 'r-/ J 


1 

2n 



1 

rc 




(19) 

The only remaining difference between the subsonic solution for the 
distribution of sources and doublets, equation (10), and the super- 
sonic solution, equation (19)^ Is the integration over surface X, 
The detailed discussion of the contribution of the surface integrals 
over X will be deferred to Appendix A. It must suffice for the 
present to remark that in all applications the integrals over X in 
equation (I9) ar© either zero or combine to give zero. It therefore 
follows that 
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fl(x,y,z) 



dxi dyi 


( 20 ) 

and the complete analogue to equation (lO) has been established throu^ 
the use 0/ the concept of finite part. 


APPLICATIOHS 

Interpretation of Previous Bes\ilts 

As a means of Indicating the various problems to vhloh equation (20) 
can be applied, three previously published results will be discussed, 
Ihese applications Include, first, the expression for the perturbation 
velocity potential of a symmetrical nonlifting airfoil (reference 6), 
second, the calciilatlon of pressure distribution over a semi— inf ini te 
wedge with leading edge swept back of the foremost Mach cone 
(reference 5) and, third, the integral equation method for determining 
the load distribution over a lifting surface of arbitrary slope 
(reference 7). 

As in the case of equation (U) for subsonic flow, let 
represent velocity potentieil $ and consider the case of a symmetrical, 
wing at zero angle of attack. 

Then = Wq and = W7 , where w„ and w, are 

az bz " ^ 

Induced vertical velocities on the upper and lower surfaceSj 
respectively. Moreover, ® ^ ® symmetrical case so 

that, since w^ — wj = 2 


n J Jr Vq 


The integral in this equation is finite at P so the finite part 
sign may be disregarded and 


« 



^ dxi 
^c 


dyi 


( 21 ) 
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This equation agrees with results given by Puckett In reference 6. 

As another example consider the solution used by B.T. Jones 
in reference 5 for a nonlifting symmetrical wing. Setting n 
equal to w in equation (20) and using the fact that w and u 
are related by the expression 


u = 



it follows that 


u = 




dz 



djci dyi 


( 22 ) 


For a wedge swept behind the forward Mach cone and having as the 
equation of its leading edge the relation Xi = myi, the 
expression for u may be written in the form 


u 



z dxi 


[(x-Xi)2 _ (y_yi)2 _ z2]3/2 


where 



Performing the Integration with respect to xi, it follows that 


u 


it Sx 


zdz 


-Yi 


- (x-^) dy. 




Uy-yi)2 + z2] 7(x - (y-yi)^ - z2 


and, after reversing the order of integration 
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u = -11 ^ / 4yi / 

n J 

y (x - ^)2-(y-yi)‘ 


- (x - ^) zdz 


C(y-7i)^+z^ y(^ - §") -(y-yi)^-z' 


= k 


2fu2. 

2n ax Jo 


In 


-y (X - - (y-yJ2 _ ^2 + k (x - ^) J 


dyi 


where k = 1 for “1 * 'TakliiS the partial 

derivative with respect to x and noting that the value of the logarithm 
at the upper limit la zero, the value of the induced velocity is 

Yi 

- 

•/(^ - - (y - yi)^ - 

and Integration yields the final result 



u = 




m 


*yrr 


In 


m2 


(x - my) (1 - m2)(x2 - jS - z2) 

y (mx — y)^ + (1 — m^) 


Zp 2u Wu 

Denoting pressure coefficient — ty — — and setting — equal 

q V V 


to 


the slope of the surface, this may te written 

(x - my) -</(!- m^)(x^~^ y^ - z^) 






In 


y (mx — y)^ + z^ (1 - m^) 


(23) 


Equation (23) gives the preasiire coefficient at any point in the 
field produced hy a wedge swept behind the Mach cone , When z is 
set equal to zero the pressure distribution over the wedge itself is 
determined and the equation corresponds exactly with eq\aation (12) 
of reference 5» 
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V/hen loading ig to be prescribed over a thin lifting aiirface, fl 
maj be assumed equal to the perturbation velocity u. A direct 
consequence of this assumption is that in equation (20) 


since, from conditions of irrotationality. 


^Uu _ ^ui 

^z bz 

By definition 

q " V ' q “ V 

and load distribution in coefficient form, , is given by the 
relation ^ 


so that 




2 u^ 2 \x^^ Au 



Equation (20) can therefore be written in the fonn 


u 1 

nr ^ A ^ 

<1 1 
•II 
1 

Jl """ 51 (it 


1 


// 


z:£!p 


z dxi dyi 


y [ (x-xi)^ - (y-yi)^ - z^]3 


(24) 


If equation (24) is transformed to the original space variables, 
the relation for u is 
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u(X,Y,Z) 

V 



2 _ p2 [(y_Yi) 2 + z2] 


Z£lp Z dJi cLYi 




3/2 


(25) 


Eqmtlon ( 25 ) is Talld for artitrary plan forms with knovm load 
di s trihut Ions . Particular examples which may he worked out with 
relative ease are the lifting surfaces carrying constant load. Once 
u is known the value of w can he determined from the integral 


and from w the ordinate z of the surface as a function of i and 
y is given hy 


where Z.e. denotes the leading edge. A discussion of trapezoidal, 
rectangular, and triangular plan forms with constant loading is 
given in reference 2 although the method of derivation is different. 

Interest in constantly loaded wings has been based primarily 
on the fact that in certain cases they can he combined to produce 
surfaces of given camber. Thus, a superposition of trapezoidal plan 
forms of variable rake, the constant loading over each trapezoid 
being a function of its rake angle, can be used to produce a flat 
plate of trapezoidal or rectangular plan fonn at an arbitrary angle 
of attack. In this case the loading as a function of rake angle is 
determined so that Induced vertical velocity is kept constant. For 
problems in conical flow a lifting element can be constructed by 
subtracting from a constantly loaded right triangle with angle of 
sweep equal to 6 the constantly loaded right triangle with sweep 
angle equal to 6 — d6. The resultant element carries a constant 
load and has a sweep angle equal to S. By summing these elements 
it is possible to find the load distribution as a function of 6 
such that certain flat lifting surfaces at angles of attack are 
formed. In reference 2, triangular wings swept back of the Mach 
cone were studied by this method for arbitrary angles of yaw. Brown 
(reference 7) has used this same lifting element to study the more 
restrictive case of the symmetrical triangular wing. 

A brief discussion of differences existing between the methods 
for producing the swept— back lifting element will shed some light 
on the various lines of attack. The approach used in reference 2 
is essentially mathematical in that a particular solution of the 
partial differential equation is used in conjunction with Green's 
theorem to satisfy the boundary conditions of the problem. The 
principal criticism of such a method is that the physical 
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Interpretation is missing. The use of equation (25), however, 
removes all such criticisms for precisely as In the case of 
Incampresslhle flow the lifting element is created by distributing 
doublets over the wing. In Brown’s solution it was necessary for 
him to determine first a line of sources by means of an integration 
aJ-ong the line and then to form the doublet line by differentiating 
along the normal to the line , The order of differentiation for 
incompressible flow is ImmaterlaJ., since the limits in the integral 
are independent of the position of the point P at X, Y, Z. Super- 
sonic flow destroys this property and it is only after the introduc- 
tion of the concept of "finite part" that the derivative of an integral 
may be written as the integral of the differential, coefficient of the 
integrand. Equation (25) thus simplifies the analysis and at the 
same time maintains the analogy with previous work. 


The usefulness of equation (20) is not at all restricted to a 
synthesis of previously known solutions. As an example of its 
generality consider its application to the problem of the rolling 
wing with leading edge swept ahead of the Mach cone. Figure 
shows the boundary conditions involved. The value of w is 
specified over the wing and, since the Mach cone is behind the 
leading edge, the value of the perturbation velocities u, v, and w 
are of course zero ahead of the leading edge. Assume for the 
moment that a symmetrical body at zero angle of attack is considered. 
It follows that if = w^ and ^ = Wj then equation (20) can 
be written in the form 


since, for reasons of symmetry, the normal gradients of w on the 
two surfaces are equal. Using now the fact that the Mach cone is 
behind the leading edge then the pressure over the upper surface is 
independent of the shape of the lower surface and eqiiation (26) 
may be applied directly to the rolling flat plate if w^ is 
detemined from the given induced velocity on either the upper or 
lower surface. This method of approach, of course, limits the 
solution to cases where the leading edge is ahead of the Mach cone. 

If the rate 'of roll is given as P 
2wu = 2PYi and equation (26) becomes 


Load Distribution for Rolling Wing 



(26) 


radians per second then 
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w 


P 


f£ 


7 

I — (27) 

T [ (X-Xi)2 - (Y-Ti)^ - I®''® 


The area t In equation (27) is that contained between the leading 
edge and the trace of the forecone on the XT— plane. Figure 5(a) 
shows the configuration for three traces corresponding to forecones 
from the points Pi, Pa, and P 3 . The region containing the 
point Pa 1 s distinguished from that containing Pi and P 3 by 
the fact that t for Pa lies ahead of the Mach cone from the apex 
and, furthermore, entirely on the rl^t of the ii ails. The 
regions corresponding to Pi and P 3 differ in the fact that ■when 
Integrating from +«o to z to find u, the upper limit of the 
Integral In the first case Is the Mach cone X^ — 3^^ — = 0 

whereas in the latter case the upper limit Is the leading-edge 


The solution must be carried out separately for each of these 
regions but only the details for the region corresponding to Pi 
need be given here since the others are similar. 

It follows that the induced velocity u at the point P is 
given as the sum of the triple integrals. 


mX - Y 


(fig. 5(b)). 


wedge Z = 







ZqYi (X + k — ) dYi 


I <yX2-p2I2 


( 28 ) 
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where A(k, Zq) is the value of Yi determined hy the intersection 
of the forecone with the leading edge on the ri^^t eind left sides 
with k equal respectively to -1 and +1. See fig. 5(a).) Ohus 


A(k,Zo) 


/ X 

-(^ m 


3^), + 7(1 + k J)2 + Zo2(^ - 3*) 



2 


After reversing the order of integration and integrating vith 
respect to it follows that 





-f 

-4/ 


A(k,Z) 


TidTiln 


-(X+k^^) - y (x+k5^)^-^^(Y-Yl)^-^V 

-(X+k^) +y (X+fc^)2-^2(Y-Yi) 2-p^2 


Moreover, since the integrand is zero at Yi = A(k,Z) the 
derivative with respect to X can he taken inside the integral and 





A(k,Z) 


Yi dY i 

'(X + ~k)2 _ 


(Y-Yi)2 - Z^ 


Integrating in this equation and camhlnlng terms it follows that 
Induced velocity u is given hy the expression 


P 2 
u = — m 


+ X 


P 2 mp^ - X 


arc sin 


mp^ + X 


. m 

It 


arc sin 


y(niX+Y)^ + (l-m^ p2)z2 
mP^ - X 


(m®p®-l)®/2 L \ 375^I^y2T7^m^P^)Z^ 



(29) 
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Setting Z = 0 In equation (29), pressure coefficient 
Ap 2 u 

— - is given ty the eapreaalon 

q V 


^ 2P r - X r mp^ - X " 

— = — -s j- arc sin ■> - • 

q ,^7 '' (m2p2 - 1) L p(mX - T) J 


mp^ + X 
(m2p2 _ 1 ) 3/2 


r mP^ + X 

arc sin 


mP^ n 


L p(mX + Y) J 


(m2p2 _ 1 ) 3/2 


,j(30) 


^rtien X>pY . 

This solution holds in "both regions containing the points Pi 
and P3 . However, in the region ahead of the Mach cone hut still 
on the wing (region corresponding to Pg) it is easy to show that 


q 


2 Pm^ 

Y(m2p2 _ 1 ) 3/2 


(mpsy - X) 


(31) 


where — <I<PY, Figure 6 shows a spanwise plot of — for 

Bi q 

m as 2 and p « 1. 

Equations (30) and (31) provide sufficient information for the 
calculation of the stahillty derivative for damping in roll, Cj . 

• Jr 

Integrati( 3 n of the load distribution yields 'the result that 


Sc, 1 

Cj =* ' — a — — 

^ S(Pb/2V-) 3P 


Load Distribution for Pitching Wing 

Another Hlmple application of equation ( 20 ) is found in the 
solution to the problem of the pitching wing. Figure 7 shows the 
boundary condition involved which 1s that the vertical induced 
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velocity "be a linear function of Xx. If Q is the rate of pitch in 
radians per second, then v = QXi on the wing. Again the solution 
is obtained only for wings which have leading edges swept ahead of 
the Mach cone. (Although solutions can be obtained for leading edges 
swept behind the Mach cone, they Involve integral equations and do 
nothing to illustrate further the direct methods of this report.) 

In the rolling wing case, ft was set equal to perturbation 
velocity w and as a result a distribution of doublets was used 
in equations (2j) and (28). As an example of the manner In which 
source— ainh distributions may be used for the same type of problem, 
equation (22) will be applied In the present case. Since the wing 
is swept ahead of the foremost Mach cone, induced effects on the 
upper and lower surface are independent and 


$ 



// 


QXi dX, dYi 


y(X-Xi)^ - (y-Yi)2 - 2^ 


( 32 ) 


Again three regions containing the points Pi, Pg, and P 3 are 
distinguished (fig. 5 ( 8 -)) and the solutions will be derived only 
for the region containing Pi. Integrating first with respect to 
Yi and then differentiating with respect to X yields 


k=-l,l 




N y-kmXi 

Xi ^ arc sin 


/(X-Xi)2 _ p2z2 

(33) 


where 


B(k,Z) = m 


X - 0^Ykm - pk y(Xmr4cT)^ + (l-m^P^) 


1 - m-'p 


2 q 2 


Considering the limit as 


(since 


q 



z — > 0 


and integrating gives: 
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''^,22,3/2 2 m I 

■ 7 / 2 Q P = ~ -mX(2^232) 


mZ(3-m^3^)-Y 


mX(2-^n^3^)+Y 



(3V) 


Formula (3^) la valid for the regions Pi and P3 of figure 5(a). 
For the region P2 the solution is: 


C~) ^ = Y - 2Xin + Xm®3^ 


(35) 


Equations (3^) and (35) provide siifficlent information for the 
calculation of the stability derivative for the damping in pitch mq, 

Yalues of 


as determined from these equations are 
plotted in figure 8 for m = 2 and 3 = 1. 


\q / \ 27 / 


Ames Aeronautical Laboratory, 
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APPENDIX A 

Discussion of Conditions on Surface X 

By definition X is the surface on which the streamlines of the 
flow first e^erlence pressure disturbances, that is, the surface 
along which the stream first becomes aware of the existence of the 
wing. Figures 2(a) and 2(b) were introduced to show the nature of 
the configurations Involved for two different plan forms. It is 
apparent that when the wing is swept ahead of the foremost Mach cone 
the wedge— like form of X is comparable to the wedge appearing in 
purely two-dimensional problems while the wing swept back of the Mach 
cone has for its siurface a conical surface and thus may be thou^t of 
as Involving a purely three-dimensional problem. 


30 


MCA TJf No. 1515 


In order to determine the value of on X it is sufficient 

to Impose the condition that the tangential component of velocity is 
continuous across X, . Such a condition represents no essential 
restriction since it is an 'Immediate consequence of continuity of 
mass flow and continuity of the tangential component of momentum 
across the surface. As a result of this condition, however, it 
follows that the 'tangential component of the pert'urhation velocity 
is zero on -the downstream surface of X since it is obviously zero 
on the upstream surface. Moreover, velocity being equal to the 
gradient of the velocity potential the perturbation— velocity potential 
must be equal to a constant on X . But an arbitrary constant can be 
added or subtracted from the velocity potential so that with no loss 
of generality the value of 4> on X can be assumed equal to zero 

and, since the conormal lies on the surface X , is also zero. 


The complete analogue to equation (10) has now been developed for 
® BO that 


ft(x,y,z) = - 


2it 







(AI) 

When Is equal to one of the perturbation— velocity 
components^ it is obvious that boundary conditions over X and t 
cannot be considered to be absolutely arbitrary since it is 
necessary to include the added restriction that the resultant 
potential ^ also satisfies the equation 


^ _ 

bz^ 

Considering first the case where the wing is swept behind the Mach 
cone, it follows that 

^2 ny 

<b= u(ii,y,z)dJCi =/ T(x,yi,z)dyi 

Vz2+y2 



w( 2 ,y,zi)dzi 
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and, after evaluating the partial derivatives of $ and substituting 
in the given differential equation, direct calculation leads to the 
conclusion that on X the following differential equations hold 


2y 



2z 



Ul 


0 


2i ^ + 2z 
bx 



Vl 


0 




Wi = 0 


where Ui, Vi, and Wi are the values of u, v, and w on X. 

The general solutions of these linear partial differential equations 
can be written as follows 


Ul 


1 

y y£+ z2 







1 

yx2+z2 


f 2 



Wi 


1 

</i2+y2 



It has been stated, however, that the tangential component of the 
total perturbation— velocity vector vanishes on X , or, in 
analytical terms 


Zui + mvx + nwx = 0 


where I, m, n are direction numbers of any tangent to X and 
therefore satisfy the relation 


Ii — my — nz = 0 


Substituting the knovn expressions for ux, vx, wx, it follows that 



or, using a different notation. 
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(A2) 


Consider now the special caae when 1*0 and m « - 22-. Iftider 
these condltlcma ^ 




+ 




so that 



Since the variables ^ and ^ are separated, the solution of this 
equation maj be written 



and 



where K is a constant. Returning now to the case irtiere I, m, and 
n are In the ratio i: 7 ; z, direct substitution into equation (A2) 

gives 







= 0 


so that 



0 


and 




MCA TN No. 1515 


33 


This equation can, however, be written In the fom 


rUf) 



y z 


from which it follows that K = 0 and Fx = Fg = F 3 = 0, All 
P®^^i^tation velocity components are seen to be zero, consequently 
equation (Al) is valid for all cases in which the wing is swept back 
of the Mach cone. 

A discussion of conditions on the surface \ will next be given 
for the case where the leading edge of the wing lies along the y-aiis 
(fig. 2(a)) and where fl represents u or w. The perturbation- 
velocity potential <b may be given by the relation 


where the plus and minus signs in the limit apply, respectively, to 
conditions above and below the xy-plane. Since * must satisfy 
the basic differential equation, an added restriction is Imposed on 
u and as a result of this condition it can be shown that 


where ui is the vaJ-ue of u on either the lower or upper surface 
of the wedge. It follows that the values of u on the two surfaces 
are 



±z 



ui “ fi(i,y) and ux » f2(i,y) 


and, since the solution Is Independent of z, and x is 
proportional to z in both cases, the final expressions are 


Ux “ f“i(y) and Ux - fa(y) 


If $ is defined as an integral involving w, the same type 

of analysis leads to the conclusion that w on the two surfaces can 
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also be eaqpreased as fxmctlons of y alone. Perturbation veloolty 
V will not be considered for this type of leading edge since the 
Inclusion of u and w covers all commonly used boundary conditions. 

It remains to substitute the results Just obtained Into equa- 
tion (19) In order to study the contribution of the Integrals over 
X . Apparently only one tenn In each Integrand need be considered 
since the conormal Is perpendicular to the y— axis and the gradient of 
9 In that direction vanishes . As a preliminary step to setting up 
the Integrals It Is convenient to Introduce a new coordinate system 
x", y”, z" which Is obtained by rotating the axial system about the 
y-axls so that the x"— and z"— axes lie respectively In the lower and 
upper wedge planes while the y”-oxls coincides with the j-axls. 

The transformation of variables Is 


When fi B u, the last two Integrads In equation (I9) may now be 
written 


x" B ^ (Xi-Zi) 



(Xl+Zi) 




1 



j •> 

dz 



di” 


Substituting for r, 


this expression becomes 
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or 


^ V—y V i2-( y— y 1 ) S-z 2 


^ oyVxa-z 2 fg(7i)clyi 

2« VyV X2-Z2 yia_(y_yj2_22 


It Is apparent that if fi(y) = — faCy) the Integrals combine to 
glTe zero so that equation (Al) may again be used In all calculations; 
moreover, the same condition applies when ft = w, Bxe aseun 5 )tlon 
that fi(y) = — fa(y) Is equivalent to postxilatlng that In all cases 
fi(y) and fs(y) are odd functions of y. In application, hovever, 
this property Is always maintained. 

It remains finally to consider the case \dien the leading edge 
of a wing Is swept ahead of the Mach cone and vrtien fl Is a 
perturbation-velocity con 5 )onent. As a means of avoiding unnecessary 
complication In treating the problem It Is possible to substitute 
first the transformation (reference 5) 


5 = — y + mx 
= — I + my 



where the leading edge has the equation y = mx, z = 0, and 
m>l. Direct calculation shows that the basic differential equar- 
tlon and the Mach cone are Invariant under the transf onnatlon and 
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that In the new ohllque coordinate system the leading ed^ lies 
along the q-ajci.s while the planes of the wedge became 

I + t = 0 


Because of the inyariant properties of the transformation, and the 
fact that the z = 0 plane is filed, equation (I 9 ) is applicable 
directly to the boundary value problem for the swept wing in the 
new coordinate system. The treatment of the integrals over X can 
therefore be developed algebraically in exactly the same manner that 
applied to the previous case, hence u and w are constant along 
the lines 


1 + ^ = 0, Tj = const 

and, again, if conditions of skew symmetry are maintained above and 
below the z = ^ = 0 plane the integration over the surfaces X 
cancel, Thua equation (20) is seen to be valid for fl equal to 
perturbation— velocity potential or perturbation velocity for all 
types of stria^t leading-edge configurations. And this is the 
complete analogue of the subsonic theory. 
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Z 



n, = cos (X i)^ - cos (7T-a) 

= cos b cos b 

Hj = COS C 7^3 - C05 C 

figure: 1- GEOMETRIC RELATIONS BETWEEN DIRECTION 

COSINES (n, ,n2,n5) of normal and direction cosines 

of conjormal to surface 5. 
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(a) rectangular plan form 

Figure 2.- Mach forecone from point P(X,Y,Z) intersecting surface t 
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Figure 2.— Concluded. 
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